Bell-type inequalities to detect true n-body non-separability 
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We analyze the structure of correlations among more than two quantum systems. We introduce 
a classification of correlations based on the concept of non-separability, which is different a priori 
from the concept of entanglement. Generalizing a result of Svetlichny [Phys. Rev. D 35 (1987) 
3066] on 3-particle correlations, we find an inequality for n-particle correlations that holds under 
the most general separability condition and that is violated by some quantum-mechanical states. 



Quantum mechanics (QM) predicts remarkable corre- 
lations between the outcomes of measurements on sub- 
systems (particles) of a composed system. This predic- 
tion is consequence of the linearity of QM, that allows 
to build superposition states that cannot be written as 
products of states of each sub-system. Such states are 
called entangled. Entanglement is at the heart of some of 
the most puzzling features of QM: the Einstein-Podolski- 
Rosen argument, the measurement problem, the paradox 
of Schrodinger cat... In the last decade, it has been no- 
ticed that entanglement is also a resource that allows to 
perform tasks that would be classically impossible. This 
new, more effective standpoint caused a renewed interest 
in the study of quantum correlations pi . 

The correlations among more than two quantum sys- 
tems have been the object of several recent studies, also 
motivated by the fact that experiments aimed to check 
such correlations are becoming feasible. Usually, the 
structure of correlation has been classified according to 
entanglement. In this Letter, we propose a complemen- 
tary classification, in terms of non-separability. Before 
entering the technicalities, it is useful to explain why the 
concepts of entanglement and of non-separability are a 
priori different concepts. We do this on the simplest 
case, that of correlations between three particles. 

The classification through entanglement presupposes 
that the system of three particles admits a quantum- 
mechanical description. Thus, any state of the system 
is described by a density matrix p. To classify a given p 
in terms of entanglement, one must consider all possible 
decompositions of the state as a mixture of pure states 
p = J2iPi l^iK^il- Then (i) If there exist a decompo- 
sition for which all are product states \ip})\ipj)\ipi) , 
then p is not entangled at all, that is, it can be prepared 
by acting on each subsystem separately. For this situa- 
tion, we use the acronym 1/1/1 QM. (ii) If all can 
be written as either \ip} 2 )\ipf) or \tpj 3 )\ipi) or \ipf 3 )\ipj), 
and at least one of the \ipf k ) is not a product state, then 
p is entangled, but there is no true three-particle entan- 
glement. We shall say that p exhibits two-particle entan- 
glement, and use the acronym 2/1QM to refer to it. (iii) 
Finally, if for any decomposition there is at least one 
that shows three-particle entanglement, then to prepare 
p one must act on the three subsystems: p exhibits true 



three-particle entanglement (acronym 3QM). 

It is difficult to establish to which class a given p be- 
longs, because in principle one should write down all the 
possible decompositions of p onto pure states. In fact, to 
date no general criterion is known. However, we know 
a sufficient criterion: there exists an operator Ai 3 such 
that: (a) if Tr(pA4 3 ) > 1, then certainly p is entan- 
gled; (b) if Tr(pM 3 ) > \/2, then certainly p exhibits 
true three-particle entanglement. The operator M 3 is 
the Bell operator that defines the so-called Mermin in- 
equality we shall come back to it later. 

The classification through non-separability (or non- 
locality) does not presuppose that the system of 
three particles admits a quantum-mechanical description. 
Rather, we have the following cases: 

(i) Each particle separately carries a script A, which 
determines the outcome of each possible measurement. 
When the experiment is repeated, the script can be differ- 
ent; the script A occurring with probability p(X). This is 
the so-called local variables (Iv) model, which we will also 
denote as 1/1 /IS. More exactly, in the Iv model the joint 
probabilities -P(ai, 02,03) that an arbitrary experiment 
A\ performed on particle 1 yields the result a\, while the 
measurement A 2 performed on particle 2 yields 02 and 
the measurement A3 on particle 3 yields 03 is given by: 

P( ai ,a 2 ,a 3 ) = f dXp(X)P 1 (a 1 \X)P 2 (a 2 \X)P 3 {a 3 \X) (1) 



where Pi(a\\X) is the probability that when particle 1 
carrying the script A is subjected to a measurement of 
A\ it yields the result a\, and similarly for P 2 and P3. 
Note that Iv is more general than 1/1/1QM. For instance, 
in Iv one can build a state such that the measurement of 
the spin of particle 1 along both directions z and x gives 
+1 with certainty, while such a state does not exist in 
QM. 

(ii) The intermediate case, first considered by 
Svetlichny ||, is an hybrid local - nonlocal model: for 
each triple of particles, we allow arbitrary (i.e. nonlo- 
cal) correlation between two of the three particles, but 
only local correlations between these two particles and 
the third one; which pair of particles is nonlocally cor- 
related may be different in each repetition of the experi- 
ment. If we define pij to be the probability that particles 
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i and j are nonlocally correlated, then in this model 

3 

■Pi,2,3(oi) 02,03) = X! Pi '3 / d^[pij(X) 
fe=l 

P i j(a i ,o i |A)J\(>4fc = Ofc|A)] ) (2) 

where fc} is an even permutation of {1,2,3}. We 
refer to this situation by the acronym 2/1 S. Note again 
that 2/ IS is more general than 2/1QM, since we don't 
require that the two correlated particles arc correlated 
according to QM. 

(iii) The last situation (3S) is the one without con- 
straints: we allow all the three particles to share an ar- 
bitrary correlation. 

It is not evident a priori whether three-particle en- 
tanglement 3QM is stronger, equivalent or weaker than 
2/ IS. The proof that 3QM is actually stronger than 2/1 S 
was given some years ago by Svetlichny Q, who found 
an inequality for three particles that holds for 2/ IS and 
is violated by QM. In this Letter, we are going to ex- 
hibit a generalized Svetlichny inequality for an arbitrary 
number of particles n, that is, an inequality that allows 
to discriminate n-particle entanglement nQM from any 
hybrid model k/(n-k)S. 

The plan of the paper is as follows. First, we intro- 
duce the family of the Mermin-Klyshko (MK) inequalities 
, that will be the main tool for this study. With this 
tool, we re-derive Svetlichny's inequality for three parti- 
cles and compare it to Mermin's. We move then to the 
case of four particles, and show that the MK inequality 
plays the role of generalized Svetlichny inequality. Fi- 
nally, we generalize our results for an arbitrary number 
of particles n. 

We consider from now onwards an experimental situ- 
ation in which two dichotomic measurements A,- and A'- 
can be performed on each particle j = 1, ...,n. The out- 
comes of these measurements are written a,j and a'j , and 
can take the values ±1. Letting Mi — a\, we can define 
recursively the MK polynomials as 

M n = 1 M„_i (a n + a'J + - M' n _ x (a n - a'J , (3) 

where M' k is obtained from Mj, by exchanging all the 
primed and non-primed a's. In particular, we have 
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E(A[A 2 ) + E(A 1 A' 2 ) - E(A[A' 2 )) , (7) 



M 2 = - (a\a 2 + a[a 2 + a\a! 2 — a\a!/) , 



(4) 



Ms — -(a 1 a 2 a 3 + aia' 2 a 3 + a[a 2 a 3 — a^a^a^) . (5) 

The recursive relation (J3J) gives, for all 1 < k < n — 1: 

M n = \ M„_ fe (M k + M' k ) + X - M' n _ k (M fc - M' k ) . (6) 

We shall interpret these ploynomials as sums of expecta- 
tion values, eg. we shall interpret M 2 as 



where E(AiA 2 ) is the expectation value of the product 
A\A 2 when A\ and A 2 are measured (note that A\ and 
A[ cannot be measured at the same time) . We call quan- 
tities such as E(AiA 2 A 3 ) correlation coefficients. We 
shall look at the values of these polynomials under QM 
and hybrid local/non-local variable models, and show 
that they give generalised Bell inequalities. 

We shall first look at hybrid local/non-local variable 
models. For technical simplicity, throughout this pa- 
per we consider only deterministic versions of the hy- 
brid variable models, which means that the script A in 
eq. (Q) and (|J) completely determines the outcome of 
the measurements (i.e the probabilities P\{X — x\\) and 
similar are either zero or one). It is known that any non- 
deterministic local variable model can be made determin- 
istic by adding additional variables In addition, we 
can also use the script A to determine which particles are 
allowed to communicate non-locally: eg. for 3 parties the 
probabilites are now given simply by 



P(ai,a 2 ,a 3 )= / d\p(\)P(a 1 ,a 2 ,a 3 \\), 



(8) 



where for each A the probabilities must factorize as some 
2/1 grouping (though not necessarily the same for dif- 
ferent A). Now, for any A, the outcomes of all products 
AiA 2 A 3 etc. are fixed, and so we can define the fixed 
quantity M„. The value of M n is just the probabilistic 
average over A of M„. Thus if we can put a bound upon 
all possible M* then we have a bound upon M n . For 
example, it can be shown that for any Iv model, M n < 1. 
This can be easily seen from (||) using a recursive argu- 
ment, noting that for any script of local variables it holds 



that either 



a„ or a r , 



In particular, M 2 < 1 



for Iv is the Clauser-Horne-Shimony-Holt inequality for 
two particles [||. On the opposite side, if we consider the 
model without constraints nS, then M n can reach the 
so-called algebraic limit M^ lg , achieved by setting at +1 
(resp. -1) all the correlations coefficients that appear in 
M n with a positive (resp. negative) sign. So for example 
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Turning to QM: since we consider dichotomic mea- 
surements, we can restrict to the case of two-dimensional 
systems (qubits) 0. In this case, the observable 
that describes the measurement Aj can be written as 
dj ■ a = a a . , with dj a unit vector and a the Pauli ma- 
trices. The equivalent of M n is the expectation value 
of the operator M. n obtained by replacing all a's by 
the corresponding a a . It is thus known that QM vi- 
olates the inequality Tr(pAi n ) < 1. More precisely, 
it is known |^| that: (I) The maximal value achiev- 
able by QM is Tr(pA^„) = 2^, reached by the 
generalized Greenberger-Horne-Zeilinger (GHZ) states 
-^(|0...0) + |l...l)); (II) If p exhibits m-particle entan- 
glement, with 1 < to < n, then Tr(p.M n ) < 2 i T J - M. In 
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other words, if we have a state of n qubits p such that 
Tr(pMn) > 2~ ~ , we know that this state exhibits at 
least (m + l)-particle entanglement. This means that 
the MK-polynomials allow a classification of correlations 
according to entanglement. But do they allow also the 
classification according to non-separability? The answer 
to this question is: yes for n even, no for n odd. As an- 
nounced, we demonstrate this statement first for n = 3, 
then for n = 4, and finally for all n. 

Three particles. Let's take the Mermin polynomial Af 3 
given in (|B|). We have already discussed the following 
bounds: Mi v = 1, Af 3 2/1QM = y/2, Af 3 3gM = M^ 9 = 2. 
We lack the bound for 2/ IS. This is easily calculated: 
consider a script in which particles 1 and 2 are corre- 
lated in the most general way, and particle 3 is uncor- 
rected with the others. Then we use @), that reads 
M 3 = \ M 2 (a 3 + a' 3 ) + \ M 2 (a 3 - a' 3 ). For any par- 
ticular script, as we said above, a 3 can only be equal to 
±a 3 . Without loss of generality, we choose 03 = <z 3 = 1, 
whence M^ 1S = maxM2. Since particles 1 and 2 can 

have the highest correlation, maxAf2 = M?/ 9 here. In 
2/15 

conclusion, M 3 = 2. Thus, for Mermin's polynomial 

Af 2/15 = Ml QM = Mf 9 = 2 : (9) 

the Mermin polynomial does not discriminate between 
the deterministic variable models 2/ IS and 3S, and the 
quantum-mechanical correlation due to three-particle en- 
tanglement. The deep reason for this behavior lies in the 
fact that M3 has only four terms: the correlations 0,10,20,3, 
a\a' 2 a 3 , aia 2 a' 3 and aia' 2 a' 3 do not appear in M3 given in 
(H). But these correlations are those that appear in M 3 ; 
thus we are lead to check the properties of the polynomial 

S3 = \ (M 3 + M 3 ) = l - (M 2 a' 3 + M' 2 a 3 ) . (10) 

For both Iv and 2/1S, the calculation goes as follows: 
we choose a 3 = a' 3 = 1, and we are left with S3" = 
\ max(M 2 + M 2 ). But M 2 + M' 2 = aia 2 + a[a 2 , which 
can take the value of 2 in both Iv and 2/1S. There- 
fore S l 3 v = S^ 1S = 1, and this implies immediately 
gi/iQM _ ^ smce 2/1QM is more general than Iv and 
is a particular case of 2/ IS. The algebraic maximum is 
obviously S 3 ' 9 = 2. We have to find S 3 ^ M . As above, 
we define an operator <S 3 by replacing the a's in the poly- 
nomial S3 with Pauli matrices. On the one hand, we 
have 

Tr(pS 3 ) = ~[Tr(pM2a a 0+Tr(pM' 2 a a3 )] < y/2 (11) 

since by Cirel'son theorem || each term of the sum is 
bounded by y/2. On the other hand, we know |l(J that 
the eigenvector associated to the maximal eigenvalue for 
such an operator is the GHZ state -i=(|000) + | 111)) . 

For some settings @, we have (GHZ\S 3 \GHZ) = V2: 



the bound can be reached, that is, S 3 ^ M = \/2. Thus 
the GHZ state generates genuine 3-party non-separability 
(non-locality). We note that, in fact, S3 is one of 
Svctlichny's two inequalities [the second inequality is 
equivalent, and is associated to ^(Af 3 — M 3 )]. 

The results for Mermin's and Svetlichny's inequalities 
for three particles are summarized in Table |. We see 
that combining Mermin's and Svetlichny's inequalities 
one can discriminate between the five models for correla- 
tions that we consider in this paper. This concludes our 
study of the case of three particles. 

Four particles. As above, we begin by considering the 
MK polynomial M4. Like M 2 , and unlike M 3 , the poly- 
nomial Mi is a linear combination of the correlation co- 
efficients of all measurements. From the general prop- 
erties of the MK inequalities M, the following bounds 
are known: M? = 1, M™ M = M 2 J 2QM = y/2, 

M 3/1QM = 2 ^ M iQM = 2v ^_ algebraic Hmit ig 

M 4 ° 9 = 4 (sixteen terms in the sum, and a factor | in 
front of all). 

Now we have to provide the bounds for 1/1/2S, 2/2S 
and 3/1 S. This last one can be calculated in the same 
way as above: using (||), we have = | M 3 (a 4 + a' A ) + 
\ -M 3 (04 — 04); we set 04 = a'4 = 1, and since we allow the 
most general correlation between the first three particles 
we have maxM 3 = Alf 9 = 2. Therefore M 4 3/1S = 2. 

One must be more careful in the calculation of 1/1/2S 
and 2/2S. This goes as follows: using (^), we have 
M 4 = \ Mi, 2 (Af 3i4 +AQ j4 ) + \ M[ 2 (M 3 ,4-M 3)4 ), where 
to avoid confusions we wrote Mij instead of M 2 , with i 
and j the labels of the particles. Now, M 3: 4 + M 3 4 = 
a 3 a 4 +a 3 a4, and Af 3 ,4 — Af 3 4 = a 3 a4 — a 3 a 4 . So if we allow 
the most general correlation between particles 3 and 4, 
these two quantities are independent and can both reach 
their algebraic limit, which is 2. Consequently for both 
1/1/2S and 2/2S we obtain Af 4 " = max(Mi, 2 + M[ 2 ), 
which is again 2 in both cases. So finally 

A/ 4 1/1/2S = Af 4 2/25 = A/ 4 3/1S = 2 < M\ QM = 2V2 : (12) 

for four particles, the MK polynomial M4 detects both 
four-particle entanglement (this was known) and four- 
particle non-separability, and is therefore the natural 
generalization of Svetlichny's inequality. 

Arbitrary number of particles. For a given number of 
particles n, we discuss only the maximal value allowed by 
QM, that is the case nQM, against any possible partition 
in two subsets of k and n — k particles respectively, with 
1 < k < n — 1, that is the case k/(n-k)S. Partitions in a 
bigger number of smaller subsets are clearly special cases 
of these bilateral partitions. We are going to prove the 
following 

Proposition: Define the generalized Svetlichny poly- 
nomial S n as 
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5,i 



ML' 



n even 
n odd 



(13) 



Then all the correlations k/(n-k)S give the same bound 
S*, and the bound that can be reached by QM is higher 
by a factor y/2: 



which is bounded by 2 _ s- because each of terms in the 
sum is bounded by that quantity. This bound is reached 
by generalized GHZ states, for suitable settings JlT| . 
Therefore S™® M = 2~r" for n odd, and we have proved 
the Proposition also for n odd. 



s nQM = ^gk 



(14) 



The tools for the demonstration are the generaliza- 
tion to all the MK polynomials of the properties of 
and M3 that we used above, namely: (I) the algebraic 
limit of M k is: M, 



alg 
k 

r kQM 



2# = M% QM V2 for k even, and 
M k WM for k odd. (Ha) For k even, M k 



M aig = 2 ^ 

and M' k are different combinations of all the correlation 



■M' k andM fc - 



M' k contain each one half 



coefficients; M k 
of the correlation coefficients, and the algebraic limit for 
both is M k 9 . (lib) For k odd, M k and M' k contain each 
one half of the correlation coefficients. These proper- 
ties are not usually given much stress, but can indeed be 
found in , or easily verified by direct inspection. 

Let's first prove the Proposition for n even. In this 
case, the QM bound is known to be S"^ M = 2~%~ . As 
in the case of four particles, to calculate we must dis- 
tinguish two cases: 



In conclusion: n-particle entanglement and n-particle 
non-separability are a priori different concepts. We have 
shown that it is possible to discriminate quantum en- 
tanglement, not only against local variable models, but 
also against all possible hybrid models k/(n-k)S allow- 
ing arbitrarily strong correlations inside each subset but 
no correlation between different subsets. Experiments 
aimed at demonstrating n-particle entanglement should 
be analyzed using the generalized Svetlichny's inequali- 
ties described in this work. Jl2j | 
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• For k and n — k even: in (||), both M k + M' k 
and M k — M' k can be maximized independently 
because of property (Ha) above; therefore, we re- 
place them by M k lg . We are left with S% = 
\M k 9 max(M„_fe + M' n _ k ), and this maximum is 
So finally 



again M^ k 



l M alg M alg k = ^ 



— k 



For k and n — k odd: in (||), M n - k and M' n 
can be optimized independently because of (lib) 
above. We have then = M®f k m&xMk = 

alg Tijalg _ n^zl 



MT\M 



k 



Thus, we have proved the Proposition for n even. 

To prove the Proposition for n odd, we must calculate 
both S* and S% QM . We begin with S*. Inserting (§) in 
the definition of S n for n odd, we find 



S n = \M n _ k M' k + l -Mi 



. M k . 



(15) 



Without loss of generality, we can suppose k odd and 
n — k even. Therefore, if we assume correlations k/(n- 
k)S, M k and M' k can both reach the algebraic limit due 
to property (lib). So S* = iM fc Q ' 9 max(M„_ fc + M' n _ k ); 
and due to property (Ha) this maximum is M"_f fc . Thus 

= 2~. Let's calculate S™® M . From the polynomial 
S n given by (15), we define the operator S n in the usual 
way. Therefore for the particular case k = 1 we have 
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are taken of the form cosOjCr^ 4 
ble choice for the settings is: oik 
a\ = ... = at, _ 
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Tr(pS n ) = ±[Tr(pM r 



10-a'J+TvipM'^O-aJ] (16) 
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2/1 QM 


2/1S 


3QM 


3S (alg.) 


M 3 


1 


V2 


2 


2 


2 


s 3 


1 


1 


1 


^/2 


2 


prod. 


1 


V2 


2 


2\/2 


4 



TABLE I. Maximal values of M 3 and £3 under different 
assumptions for the nature of the correlations (see text). The 
last line is the product of the two values. 
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